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A statistical ansar. ble of large eddy sinuktions is run sir. ufcaneousV fcr the a e 
ow . The inform ation provided by the di ersit largeroale velocity eHs is usd xi an 
aasEm ble-averaged version of the dynam ic m odel T his produces localm odelparam eters 
that only depend on the statistical properties of the ow . An mportant property of 
the ensat. bl&averaged dynam it procedure is that it does not require any spatial aver- 
aging and can thus be used in folV iihccn ogeneous ows. A iso, the arsan bde of LE S 's 
provides statistics of the largctscafo vebcity that can be used for building new m odels 
for the subgrid -scade stress terror. The arsan bleaveraged dynam ic procedure has bear 
in plan aited w ith various m odeis for three cw s: decayitg isotropic turbulsnce, forced 
irotropi: turbuiaroe, and the tin e-developing plane wake. Itis found that the reailts are 
aln ost indepardart of the num ber of LE S 's in the statistical arsan bde provided that the 
ensan bfo contains at least 16 rtalisations. 


1 . Introduction 

T he num ber of degrees of feeedcm needed to draracterize a velocity edri u, that corre- 
cponds to a turbuknt ow is know n to increase as Re:' 1 * ( Re. is the R eynodds num ber) n 
three din arsional turbulent systems. D irectnum ea±alf ubtions (DNS) oftheNavier{ 
Stokes aquations governing the evolution of such system s are thus din ited to m oderateV 
an all R eynoids numbers. There is thus an interest in developing techniques m which 
only a faction of the totalnumber of degreesof faedcm tactually sm u dated. Am ong 
these techniques, Large Eddy Sinulation (LES) and Reynolds Averaged Navier{ Stokes 
(RANS) smuhtion have attracted m uch interest in the past few decades. InLES, the 
num ber of degrees of freedom is reduoed by using a spatial Iterirg: 

Ui (x) = J dyC (x - y)u, (y) , 

where (7 is the Iter kernel and tx, is the LE S eld . InR AN S , an arsan ble averaging is 
U 33 d to de ne the RANS eld U i : 


Ui = M 


( 12 ) 
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Inboth cases, the equations for 77, or for Ui contain an unknown stress term that requires 
m odelling . T he purpose of the approach developed here is to can hire cona^ts fern the 
two m ethods to produce a statistical version of LES . 

T he present approach is m otivated by the foot that, in both LE S and RAN S , m odels 
for the degrees of heed cm that have bear elm hated are inspired hem statistical theories 
of turbulence. I tis thus in plicitiy assum ed that the Jtering and ansan ble averaging can 
both be regarded as pro jective operations that associate a num ber of di erent velocity 
realisations w ith a single LES orRANS eld. There is, how ever, an in portant di ermoe 
between on sen ble averaging and ^>atial Tiering. The ensan bh averaging opera tion re- 
duces the num ber of degrees of freedom by so m uch that aln ostno useful inform a tion on 
the uctuations Suj = it* — U t can be deduced form the knowledge of [/, alone. On the 
other hand, in LES the statistics of the unresolved scales u ■ = m - 77, m ust be closeV 
related to the statistics of the resolved scales u l since there is no char scale separation 
between then . H ence, knew ledge of the statistical properties of the LES elds 77, should 
be he^fol in developing LES m odels. The advantage of studying a statistical organ - 
ble of LE S 's is the ability to extract statistical inform ation for building m odels for the 
unresolved scales. T his w ill be explored in Section 3 . 

T he de nition of equivalent and independent LE S elds is not necessarily obvious and 
dnould probably depend on the m otivation for the sin u la tion . W e propose in Section 2 
son e conditions under which two simulations of a turbulent aw will be supposed to be 
independent and equivalent. In Section 3 , w e w ill draw that the knew ledge of an ens=m ble 
of L E S 's yields a good fram ew ork for developing a local version of the dynam ic procedure 
in which m odelparam eters are can puted usiig statistical quantities. The application of 
this approach to isotropic turbulence is presented in Section 4 . Application to the wake 
aw is presented in Section 5. In this last case, it is drawn that the knowledge of an 
ensen ble of realisations can be used to develop new m odels that explicitly incorporate 
averaged quantities m ade available through the an sen ble. 


2. Statistical ensemble of LE S 's 

T he equation for large eddy sin ulation (L E S ) is obtained by app lying a ^atial Iter to 
the N avier{ Stokes equations. T he LE S equation thus describes the evolition of a Itered 
velocity eld 77, , which explicitly depends on the an all scales through the subgrid scale 
stress Tj j = TTTuJ — 77, Hij : 

<)tUi + OjUjTti ~ ^,,V 2 77 i - djT i} . (21) 
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For sin plcity, we only consider foccm pressible cws, in which p, the presaire divided by 
the density, is detenu hed by the in can possibility condition . The unknown tensor r i: 
appears in the equation for the forge-scale velocity Tii but it depends on the an all-scale 
velocity e3d . T he purpose of this study is to explore the advantages of sin uitaneousV 
running several stetistically equivalent and independent LES's for the sam e cw . In 
practise, w e thus replace the equation (2d) by the folk*/ fog set of equations for R forge- 
scale velocity elds u\ : 

dtUj + djUjU 1 ' - - djj f + v { |V 2 u^ - dj r£ , (2 2) 

where r - 1, . 

It is worth m entionfog that the use of an ansan ble of LES's is not per semuch m ore 
expensive than the use of a single realisation .To dnc m this, let us consider a stationary 
LES and denote by t t the tin e of the transient period between the beginning of the 
sin ufotion and the tin e at which the turbulence beoom es folV developed. Let us also 
denote by t s the tin e {beyond t t ) required to converge the statistics. Than, the C PU tine 
required for obtaining converged statistics w ith a single LES is t t + t s . W ith an ansan ble 
of realisations, statistics are accum ufoted over both the ansan ble and tin e. T hus, for 
equivalent sam pie, the ansem ble only needs to be advanced in tin e by the am cunt t s /R - 
T he total C PU cost for the ansan ble is thus R {tt + t s / R) , which am ounts in an overhead 
of (R - 1 )t t over a single realisation. If the ratio between the transient phase and the 
tin e needed to converge statistics is sn alb than the additional cost w ill be m oderate. 
In the exam pfos treated below, this additional cost is totall/ negligible. M oreover, if the 
LE S is not stationary and if there is no direction of hem oganeity, the ansan ble-averaged 
approach is presum ably the only way to obtain statistics. 

21. S tatis tkally oquivafont and independent LES 's 

The know fodge of an ansan ble of LES's can on V be usefol if the LES elds are all 

independent. Yet all these elds have to correspond to the sam e experin antal situation 
if son e m eanfogfol statistics are to be extracted from the ansan ble. W e therefore m ust 
de ne what will be considered statistical!/ equivalent but independent LES elds. Al- 
though a proof of existence and uniqueness of seditions for the N avper { S tokes equations 
isnot yet avaifobfo, frem a practicalpofot of view a cm described by the N av per (Stokes 
equations or by an LE S equation is assum ed to be folly de ned by the know ledge of 

1 . T he dan afo V in which the cm is considered . 

2. The cond itio ns on the boundary dT) of this donah {d'D^ t) = b t it) where the 


functions bi {t ) are given . 
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3 . T he initial conditions u, (x, 0) = u { - (x) Vx 6 V . 

H cm ever, in a sin ubttion of a turbulent cm only tine dan ain and tine boundary conditions 
are rigorously xed . Indeed, because of tine lack of sensitivity to initial conditions in 
turbulence, di erent sin ulations w ith di erent initial conditions dnaring son e properties 
are considered to characterize the sam e cm . T hus, the requiren ant that the initial 
conditions are kncwn is scr\ e^hat rebxed for turbulent avs and the point (3) is thus 
replaced by a weaker constraint: 

3 The initial condition Ui (x, 0) = uj 1 {x ;w \ ) is generated using randan numbers W( 
and satis es a certain num ber of constraints: P s [u l - ] = p s , s = 1 , . . . S . 

For exam pie, in ham ogeneous turbulence, the m ost in portant constraint w ill be on the 
qpectzum E {k) of : 

k 2 j d (k, 0)| 2 = E(k) (2 3) 

where ti { - is the Fourier transform of u* 1 and is the solid angle in the wavenumber 
= k 2 d dk . For channel cm , one could in pose the plane-averaged pro le of both 
the velocity U iy) and the R eynolds stress R l3 (y) : 


<«?),., - u t-M.i 

(2.4) 

( («" - V (y)(5,,i) 2 >i,z = X, (>/) 

(2 5) 


where x , y and z are respective^ the longitudinal, the wall norm al, and the transverse 
directions and {• * -} Xi2 represents the average in planes parallel to the wall W e w ill not 
discuss in detail the m in in al constraints that m ust be in posed on the initial conditions 
in order to have a reasonable sin ulation. In feet, this m in in al set of constraints w ill 
probably depend on the type of cw as well as on the quantities that are m ensured in 
the sin ulation . H ere we only suppose that these constraints do exist in order to give a 
precise de nition of equivalent LES 's: 

i) Two LES 's are statistically equivalent if the don ain of the ow and the boundary 
conditions are the am e and if the initial conditions satisfy the sam e set of constraints. 

Carrying an ensemble of equivalent LES's can be con putationalfy e ective onfy if 
the di erent m en bers in the set of LES's are independent. Here again, the de nition 
of independent LES 's m ight depend on the cm as well as on the quantities that are 
m easured in the sin ulation . A long the sam e line as for the de nition of statistically 
equivalent LES 's, we propose the follcw ing de nition : 

ii) Two LES 'sare statistically independent if their initial conditions are generated 
w ith un correlated random num bers W [ . 
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W e nan ark that for a stationary ow , such equivalent and independent iiitkl con- 
ditions can be obtained by runniig a single LES and recording several velocity elds 
s^arated by at hast one large-eddy turnover tin e when turbulence is folly de/eloped . 

2 2. V niversalm odelparam eter in LES 

C Irssical closure strategies in LES am ount to m odeHkng t t j ii term s of the resolved 
velocity eld : 

r i:i = C m£ 61/ , 3, ( 2 - 6 ) 

where is the Iterw iith . T he tensor is supposed to characterize the dependence on 

both the iter width and the speci c realisation of the large-scale ow tip 0 n the other 
hand, we w iUassum e that the param eter C droids onV on the type of cw and on the 
Iter drape and dnoull not depend on any particular realisation of the large-scale velocity 
eid.f In the following, we will refer to this assum ed property as the universality of the 
model param eters in LES: For a giver geometry and for a given Reynolds number, the 
m odelparam eters should be the s am e in all equivalent LE S 's. T his concept of universality 
does not in pV that the m odel param eters are constant in ^aoe and tin e. C learly, 
C - C (x, t ) can be a eld quantity that needs to be adapted both in space and tin e to 
the local conditions of the cw . H cwever, ir our approach the variations of C are not 
supposed to take their origin in possible uctuations in the large-scale cw - R ather, ( 
is expected to depend on\r on the averaged properties of the ow , and in that sense it 
shares m any properties w ith RAN S quantiti e s. 

T he assum ption that the m odel param eters are universal has a direct in uance on 
the formulation of m odels in an ensan ble of statisticnlV equivalent LES's. In the equa- 
tions (2 2), these m odels diould have the follcw hg structure: 

r[j = C m[j ], (2.7) 

where C is now independent of the realisation index r . 

I tm ust be noted, however, that the developm ant of the dynam ic procedure in som e 
ways chalfenges this viewpoint. In the dynam ic procedure, inform ation from the an all 
scate of is used for estin ating the m odel param eters. This procedure is kncwn to 
produce highly uctuating m odelparam eters. Such a property is son etin es regarded as 
a proof of the capability of the dynam ic procedure to produce m odelparam eters that 
account for the local conditions of the cw . H cw ever, these uctuations in ( are reqoon- 
sdbfe for instabilities, and scm e averaging procedures are used to avoid this di culty. 

f Of course, more sophisticated models with more than one term have also been proposed, 
but, the specific roles of the model parameters and of the model tensors remain the same 
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W e propose in Section 3 an approach that reconciles the dynam ic procedure w ith the 
concept of a universal m odelparam eter. In this sense, it is fairly di erent fcm other 
procedures in which the oonoqpt of a universal param eter has not bean adopted, such as 
the local dynam ic procedure developed by Ghosalet al (1995} , the Lagranghn dynam ic 
procedure proposed by M eneveau et al (1996), or the tin e lagging procedure proposed 
by P ion elli and L ii (1995). 

2 3 . N ew m odelling concepts 

The knew ledge of an aisanble of LES elds opens new possibilities in the m odelling of 
the r[j . Indeed, it is now conceivable to introduce an explicit dependence on ensanble- 
averaged quantities into the m odels for t/J- . 

2 3 1 . M odel has ®1 on the uctuating s tain tans or 

T he rstm odelw e proper is based on the uctuating part of the rate-of-strain tansor: 

t[j = -2v, (s\ 3 - (S i3 )) = - 2 v r 6S'j , ( 2 . 8 ) 

where u e is the eddy viscosity. T his form ulation has seen e nice properties. T he averaged 
total disspation is given by 

€= <MV‘ s 'y> + •'»(■%■% >, (23) 

and consequently the turbulent disspation originates only from the uctuating part of 
the strain tansor. Them ean part contributes only to them o lecu hr disspation .T his prop- 
erty ensures that the m odelw ill not produce disspation in a hm inar region . In addition , 
while this m odel is disspative on average (provided the eddy viscosity is positive) , indi- 
vidual realisations can have negative dis sp ation, thus representing the inverse transfers 
of energy &m the sn all unresolved scales to the huge ones (badcscatter) (Leith 1990; 
M asen & T hem sen 1992 ; C aratiet al 1995a } . I tis generally believed that backscatter 
originates horn uctuation phonon ana in the subgrid scales, and np>resantation of this 
e ect through uctuations in the strain tenser is thus very reasonable. 

R esults using this m odel for the wake cw are given in sections belcw . I thas already 
bean used in the channel cw , where the plane of hem ogeneity is used to com pute the 
average (Schum arm 1975} . H aw ever, the arsen ble of LES 'sallows the use of such m odels 
even in fully inhem oganeous cw s. O f course, m any other m odels m ight be considered 
along the sam e lines, and the uctuating strain rate is not the only quantity that could 
enter the m odel In this paper, we w ill restrict our investigations to the m odel (2 .8} in 
the study of the wake cw . H ewever, we m antion hereafter another possible use of the 
know tedge of an ansen ble of LES 's in the case of anisotropic cws. 
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2 3 2. Anisotropic m odel 

Anisotropic e ects are aln ostuniversalV observed in turbulence. Hew ever, anisotropy 
usual!/ orginabes firm can plex interactions betw een cw direction , solid boundaries and 
external constraints like pressure gradient or global rotation . Itis thus quite di cult to 
predict a priori the m ain consequences of this annsstropy. In the context of statistical 
averaged LES, we have access at any instant to m ean quantities that w ill display the 
anisotropic structure of the turbularae even for folly inhom ogemeous cws.A m od el that 
would direct!/ take advantage of the ensan ble of LE S 's could be: 


Tij ~ LHikljiSh > (210} 

where the factor /x plays the role of an eddy viscosity but through an anisotropic relation 
betw een the subgrid scale stress and the strain tensor . T he tensor 7g should be a m ©asxre 
of the anisotropy. 1 1 could be constructed firm the velocity uctuations: 

Suj) 


lij - 


( 211 ) 


<<K Su l) ' 

This m odel reduces to the chssbal eddy viscosity m odel for isotropic turbulence ( 7 i 3 = 
6 t j ) . The -srjn of the disspation droids onV on the sis 11 offi shoe the product of 77 
and the strain tensor is given by 




( 212 ) 


M oreover, if there is no turbulence in one direction (Su a = C ) , the m odelhas the property 
that the con p on arts r 1(l = r aj = 0. This is an expected property that is m jssed by 
isotropic eddy viscosity m odels. 


3 . C oupling the dynam ic procedure and the ensemble of LE S 's 
3 1. C hssial dynam ic procedures 

T he dynam b procedure is based on an hartity (G erm ano 1992 ) that relates the unknow n 
stress garerated by di erort Iters: 

Lij + Tij — Tij = 0 , ( 31 ) 

vjhereTtj = izfSj-S, % is the subgrid scale stress generated by the successive app hcation 
to the velocity eh of two Iters that are respectively denoted by - and by ~.The Leonard 
tensor is given by L {j = uT&j ~ ^ i §? - Itdepends onV on u x so that it does not require 
any m odellhg. This hortity (3 1) is of course onV valid for the exact and unknewn 
subgrh scab stresses. W her m odels are used, t 13 ~ Cmjj P/ ] and T tJ « CmJj \Pi L the 
di ereroe E rj = L tJ + CrnJj - CrriL betw ear the tight hand side and the left: hand side 
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of [31) can be considered as a m easuxe of the perform anoe of the m odeL T he dynam ic 
procedure uses this m easure xi order to prescribe the m odelparam eter C by m in in izing 
E t j . W her a hem oganeous direction exists in the p rob Jen , the estfcu atbon for C is gi/ai 
by (Germ ano etal 1991; Lily 1992;Ghosalet al 1995): 


('i ^ 

{ Mij M ij } fi 


(3 2) 


where A Uj = mjj - mj^ and the average {■ • -)h is supposed to be taker over the ho- 
rn ogaieous direction (s) . 0 bviousiy this approach is restricted to qpecialgeom etries w ith 
hem ogaieous direction (s) - C cm piex geem etries require an alternative treatrn ait in w hich 
a local de notion of the param eter C can be proposed . This is the case in the local dy- 
nam ic procedure developed by G hosal etal (1995) as well as in the Lagrang ian dynam ic 
procedure proposed by M aieveau et al (1996) . In both cases, the m odelparam eter is 
directly related to the large-scale eld Ui through the terser L u and M ?J . I twill thus 
vary from one realisation to another, evai if the underlying LE S 's are supposed to be 
equivalent. A s already m eitioned in Section 2 2 , the dynam ic procedure thus produces 
m odelparam eters that are not universal In the early stages of its developm ait, the fact 
that the m odelparam eters are directly related to the ^>eci c realisation of the cw was 
considered advantageous because this allowed the m odelto bem ore adaptative. H cw eter, 
this property proved to be problem a tic because it generates highly variable m odelpa- 
ram eters that cause num erical instabilities. Scm e of these practical problem s have beai 
resolved in the aforan aitionod local and Lagrang ian versions of the dynam ic procedure. 


3 2 . Enson ble-aveaged dynam ic procedure 


T he eisen ble-averaged dynam ic procedure (EADP) w e propose here is con cep tualV very 
close to the volum e-averaged or plane-averaged versions (3 2) of the dynam ic procedure. 
T he on V di ercnce ccm es from the nature of the average, which is new an aisan ble aver- 
age over the set of LE S 's. C onsddering that R LE S 's (2 2) are com puted sin ultaneously 
the m odelparam eter is now givai by: 


r ,_ (LiiMn) 

' ~ (MijMtj) ’ 


(3 3) 


where (• • ■} now r^reseits the aisen ble average. T he expression (3 3) is only valid if the 
param eter C is slew V dqpaidait on qpace and can be takoi out of the test Iter . Such 
an assum prion is not very restrictive, however, since the aisan ble averaging is likely to 
an ooth out the rapid variations in the eld. In the next Section, it will be seal that 
the m odel coe c:eit does indeed becom e an oother and an oother as the aisan ble size 


i 



Statistkaloisan ble of forge eddy sin ufotions 9 

is increased . T he form ufotfon (3 3 ) guarantees that the m odelparam eters are universal 
shoe they depend on V on the statistical properties of the large-scale velocity elds. 

Ina sense, the ensemble ofLES's corresponds to an arti. cial direction of hem ogeneity, 
whih alvays exists independent of the ccm p laxity of the ow . In the unexpected cases 
xi whrh the m odel ooe ciant ran ahs signi cantly variable in ^>aoe for large ansa n ble 
sizes, the E A D P could be coupled w ith other approaches like the localdynam ic procedure 
developed by G hosnlet al (1995) , the Lagrangfon dynam ic procedure proposed by M en- 
eveau et al (1996) , or the tin e Egging procedure proposed by Piom elliand Liu (1995) . 
The coupling of the EADP with any of these m ethods would lead to a negligible cost 
since it would be used only once for the whole an sen ble. 

In son e cases, the feet that the m odel param eter cannot be adapted to the speci c 
realisation of the large scale cw m ight be considered as a drawback of the EADP. In 
p artic u lar , a greater adaptability m ight be desirable for very in term itbant cm s w ith, 
for instance, fo calize d turbulent spots appearing inside a lam inar sea (H anningsonet al 
1987) . Indeed, when the turbulent spots appear random a statistically hem ogeneous 
dom ain , the m odel param eter predicted by the EADP is quasi-constant and is w eakly 
a ected by the tmfoufont spots. In feet, the EAD P in plicitV assum es that the m odel 
tensor m y Ep dioufo, alone, take care of the turbulence activity. W e realize however 
that, because the perfect m odel is not available, the assum ption of a universal param eter 
m ight be som etin es inappropriate, depending on both the nature of the ow and the 
m odel adopted for ifiij P/ ]. However, it must be noted that the sam e di culty would 
be encountered for m odefo that use param eters that are chosen a priori as w ell as for 
dynam ic m odels that are based on vofom e or plane averaging. M oreover, in the test 
cases presented in the follow ing sections, the predictions of the EADP are encouraging. 
M ore fo calized approaches (G hosal et al 1995; M eieveau et al 1996; P xm elli & Liu 
1995) for whrh the assum ption of a univerralm odelparam eter is not adopted would 
probably respond m ore strongly to in term ittant cws. Itmustbe noted, however, that 
these generalized dynam ic procedures are usually can bined w ith m odels that are basad 
on statistical approaches. In these cases, use of the focalized dynam i c m odel in a sense 
contradicts the underlying statistical assum ptionsused to build the eddy viscosity m odel 

F in ally, w e ran ark that the m odelparam eter is the only coupling betw een the di erant 
LES's. The di erance between DNS, LES, and an sisanble of LES s coupled through 
the E A D P is illustrated in gure 1 . A s a consequence, the E A D P is perfectly suited 
for distributed processing on parallel ccm puters. T he m ost natural in plan antation of 
this procedure am ounts to running each m on ber of the ensan ble of LE S on a separate 
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a DNS 



b LES 



c : Ensemble LES 

FIGURE 1. Tile differences between DNS, LES, and ensemble of LES’s using the EADP are 
illustrated. In DNS (top), only the right hand side of the Navier-Stokes equations is needed 
for advancing the velocity field in time. In traditional LES (middle), an additional modelling 
term is needed. In the EADP (bottom), one substep, common for all the LES’s, is added for 
computing the model parameters used in each of the simulations. This is the only point where 
information is required from the other fields. 


node. C cm m unication between the di erect processes is lin ited to the cam putatbn of 
the m odel coe ciait. O therw ise each eld u\ is advanced in tin e ind^aadantiy of the 
others. T his property foould guarantee very good scalability if large ensan ble sizes are 
explored. 


4. Tests in isotropic turbulence 

4 1 . D ecaying turbulaice 


The EAD P described in the previous section was tested in decaying isotropic turbulence 
for 32 :i LE S 's. T he taisor mh was chosen to correspond to the Sm agorfodcy m odel 


■■-2C\S T \Si 


(41) 


A serhs of num erical experin aits has determ hed (C arati et al 1996) how large the 
aisan ble of sin ultaneous LES 'sm ust be ( i.e . haw large/? should be) . T he criteria used 
to deteam foe the m foin al size of the aisan ble w ere focused on 
1 . T he qpatal variability of C . 

2 . T he percaitage of negative C . 
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Figure 2. Tv pica] oiie-dinieusional profiles (left) and probability distribution f\inction (light) 
of C in decaying isotropic turbulence for R = 1 (long-dashed line), R - 4 (short-dashed line) 
and R = Hi (solid line). The averaged value of (C) ~ 0.02 before clipping is almost independent 
of the ensemble size. The expected “smoothing effect” of the ensemble averaging is reflected by 
a rapidly decreasing deviation <r 2 c = <(C - (C» 2 > with R. For instance, a c ~ 0.3 for R = 1, 
o r ps 0.00 for R = 4, and ac ~ 0.02 for R = 16. 


3 . C on parson w ith the volim e-averaged dynam ic m odel 
4 . C cm parison w ith direct num erical sin ulatbns. 

The resuits are quite encouraging . 1 1. appears that with only lb sin ultaneous LES's, 
the oosan ble-averaged dynam ic m odel perform s as well as the volim c averaged m odel 
The ^atial variability of C decreases drastically whaa R h creases. T his is also re ecbed 
in the probability distcfoution fijnction (PDF) of C {^e gure 2). Seme quantitative 
m easurm acts of the spatial variability as a Sanction of the ensemble size are given h 
Tabfe 1. In particular, the faction of negative C before cipping drops fan 41% for 
/? = 1 to 15% for R = 16.Haaoe, the faction of points for which C has to be clpped is 
still sxjni cant evm. fori? = 16. However, the consequences of this clppfog are less and 
apni cant because the clpped valies of C have sn aller m agnitudes for hcreashg 
arsan bte sizes . For instance, the ratio between the averaged values of C 1 before and after 
the clipping is only 0 J.9 for R = 1 while it reaches 0 31 for R - 16 . H aioe, the e ect of 
cipping on physical quantities like the energy disspation beccm es sn all for R > 16 . 

T he ccm parison between a 512’* DNS and the dynam ic m odel shew s good agrean ait 
both for the total resolved aiergy and for the ^ectra. The aisan ble-averaged results 
for R = 16 are hd istingu idiab Je firm the volim e-averaged values and on V the EA D P 
results are plotted on gure 3 . A n aisan ble of hdqpardait volim e averaged LE S s w as 
run to albw ccm parison of both the m eans and the standard deviations. 

4 2 . Forced turbulence 

W e have run an ensan ble of 16 32 :{ forced turbulence LE S 's w ith zero m olecular viscosity . 
F igure 4 drew s that the m ean resolved energy and the standard deviation evolve sin i- 




12 


D . C arati, M . M . R ogers , and A . A . Wr ay 


R 

Before 

clipping 

After 

clipping 


(C) 

fTC 

(C) 

oc 

1 

0.018 

0.29 

0.089 

0.19 

2 

0.020 

0.12 

0.048 

0.081 

4 

0.020 

0.057 

0.031 

0.040 

8 

0.019 

0.031 

0.024 

0.024 

16 

0.018! 

0.020 

0.020 

0.017 

32 

0.018 

0.013 

0.019 

0.012 


Table 1. Average and standard deviation of the model coefficient (before and after clipping) 

versus the ensemble size. 



Figure 3. Comparison of the energy decay between the truncated DNS (solid line) and the 
averaged energy predicted by the set of LES’s using EADP (dashed line). The dotted lines 
correspond the averaged energy ± one standard deviation as predicted by the set of LES’s using 
EADP 


~hrV for both the volim e- and on sen ble-averaged m odels. T his shew s that the coupling 
induced by the can putation of the m odelparam eter through the EADP approach does 
not introduce sp urfous corrections between the di eraitm an bets of the eisan ble. T he 
standard deviations ran ain sin ilar in the two approaches, indicating that the LE S elds 
in the EADP rem ain nearV as independent as those in the ansanble of independent 
volim e-averaged sin ulations. It is also interesting to can pare the con pen sated energy 
oestrum E (k) = E {k ) k r ^ :i ( , where E {k) is the energy ^ectrum and c is the dis- 
spatbn rate. 0 f course, with 32* LES, we do not expect to observe a well developed 
inertial range or to obtain a very good estin ate of the K oin ogorov constant. H cwe/er, 
exam inatbn of the results in gure 4 indicates that the observed \K oin ogorov constant" 
is reasonable. 


5 . Tests in w ake ow s 

The ow considered here is a tin e-evolbng plane wake for which data fan both 
direct num erfoal sin ulations (M osar & R ogers 1994 ; M osar et al 1997 ) and large eddy 
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Figure 4. Resolved energy (left) and compensated energy spectrum (right) in forced isotropic 
turbulence: ensemble- (solid boxes) versus volume-averaged (open boxes) dynamic procedure. 
Dotted lines and dashed lines correspond to the mean ± one standard deviation in the FA DP 
and in the volume- aver aged dynamic model, respectively. 


sin ulations (Ghosal & Rogers 1997} are available. This cw is both statistical!/ non- 
stationary and xihcrn ogmeous and dnould thus be a m ore dan anding test of the E A D P 
than the hom oganeous cm s studied in the previous section . 

T he pseudoqo ectral direct num ericalsin uhtion of the plane wake considered here has 
bean described in detail by M osar & R ogers (1994) and M oser, Rogers & Ewing (1997). 
T he qDatialdep>endaoce of the independant variables is r^resai ted in the periodic stream 
w ise and spanw ise directbns by Fourer base Sanctions and the cross-stream dependence 
xepresanted by a class of Jacobi polyncm ials on a m apped in nite dcm ain . Up to 
512 x 195 x 128 m odes are required to accurate!/ resolve the turbulence. T he R eynolds 
num ber based on the integrated m ass ux de cit, 

fi= - f {U (y) - U^)dy, 

J — no 

& Re = fi/u = 2000. In a tin revolving plane wake, the integrated m ass ux de cit is 
constant. 

LES's of the sam e ow using the dynam ic procedure, w ith a IteredDNS eld as 
an initial condition, have been reported by Ghosal & Rogers (1997) . The emulations 
were psEndospectral like the DNS, but the spatial depandance of the vorticity in the 
imhom oganeous cross-stream direction is represented in term s of Fourier m odes on a nite 
dcm ain . T he appropriate non-periodic velocity eld is than calculated using the m ethod 
of C orral & Jin enez (1995) . The num ber of m odes used in the LE S 's w as 64 x 48 x 16 
and the sam e num ber of m odes and sam e num ericalm ethod have bean adapted (C arati 
& R ogers 1998) for the EAD P LES's exam ined here. Thus each LES requires up to 260 
tm es fewer m odes ocm pared to the DNS. 
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5 1 . The sufopdd -scale m odels 

In the present study, we have investigated three di erant m odels, all based on the eddy 
visoosity concept. T he rst one is the Sm agorindcy m odel introduced in the previous 
section (41). I nth ism odel, the inertialrange scaling for the eddy visoosity Vt ~ V 1//,J 

is expressed in term s of the resolved strain ^ate tensor by using the approx in atdon for the 
disspation rate e ~ . This approxin atbn is required in traditional LES because 

a separate equation for the disspation rate is not usualV ocm puted . Howwer, in LE S 
based on the dynam ic procedure, the din onsional product C f ~ C can be predicted 
directly- This has motivated (W ong & Lilly 1994; Caratiet al 1995b) an alternative 
m odel based on the inertialrange scaling such as 

M odel A : rfj ~ -2(7, 4/1 S’ u . (5 2) 


F in ally, w e have considered a third m odel for w hih the tensor r£ is given by the expres- 
sion (2.8) 


M odel B : t[j as -2 C € 4/4 ^9^ - (5^-)) , 


(53) 


where the brackets indicate ensam ble-averaging over all realisations. T he possible advan- 
tages of this last m odel have been discussed in section 2 3 . 

I nail three m odels, the sign of C (or ofC e ) w ill also determ ine the sign of the subgrid- 
scale disspation , since a negative C corresponds to a negative eddy viscosity. In order to 
avoid num erical instabilities, the m odelparam eter m ust than be set equal to a m in in al 
positive valie (cipping procedure, see Ghosalet al (1995)) at points where the total 
viscosity (eddy plus m olecular) is negative. For the Sm agorinsky m odel the stability 
condition 

— -2 / — r — r \ 

C (2 S kl S kl J +i/o> 0 (5.4) 


dpends on the realisation . T his is an undesirable property since ( J is supposed to be a 
universal cw characteristic for all m an bers of the ansenble. An alternative formulation 
in which C is indeed the sam e for all realisations results firm the follow ing stability 
condition 


C 


— 2 



lA\ > 0 . 


( 55 ) 


In the lin it of an in nitenum ber of reals. tions, them axin um of the resolved strain -rate 
tenser am plitude would be aln ost unbounded . H enoe, for the Sm agorinsky m odel it is 
reasonable to sin ply in pose C > 0 . Form odelA , hew ever, the situation is di erant. T he 
stability condition is naturally the sam e in each realisation 


C f 


—4/:* 


> 0 , 


I 


(5.6) 
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For sin p Icily, the sam e condition has bean used for m odelB . 

5 2. T he initial conditions 

I n practice, initial conditions for LE S can be bu lit either by IteringaDNS {whan it is 
available) or by generating a random velocity eld satisfying son e constraints {as dis- 
cussed in Section 2 1) . For the EAD P , we have to generate/? equivalent but independent 
elds. In the r^.«^ of isotropi: turbulence, the only constraint that had to be satis ed by 
the initial eld was the energy qpectrum . W e thus have used Rogallo's approach (1981) 
to build R in itial conditions w ith the sam e spectrum and independent phases. For the 
tin e-evol/ing plane wake, random initial conditions could be generated following the 
sam e approach as the one used for initialising the DNS. However, for the plane w ake, a 
large num ber of quan tities are m easurad and any num ber of then m ight be considered as 
constraints that need to be m airntained by all realisations {eg. pro les of m ean velocity, 
turbulant kinetic energy, enstrophy, etc.) . The m ain purpose of the present study is to 
test the EADP rather than evaluating an initialisation procedure. For that reason we 
have u sed a sin pie trick to generate R statistical^ identical initial elds. 0 ur procedure 
fo based on the foot that the observed quantities are cam puted through plane averages 
and are thus invariant under the change 


ui Ci’, ty, 2, Ui ) — * Ui Ct + V, z + $ z , in ) - 


(5.7) 


Thus by using R values of R ini tia l velocity elds are produced that clearfy 

satisfy the requiren ant that the LES realisations be statistical^ equivalent. H cwever, 
this procedure does not produce statistically independent initial conditions, even w ith 


random choices for because the two elds are identical and sin ply knitted in 

^paoe. W ithout the subgrid m odel term s, all the statistics would ran ain identical for 
^11 tin es. H ow ever, the m odel term s w ill have the desirable e ect of de-correhting the 
di erent m on bers of the eisen bis . T his results because the universalm odel term s act at 
the sam e C x - z) location in all the realisations, not at the sam e relative position in the 
dni&ed ows. Exam pies of this de-correlation are given in gure 5, in which the reduced 
m ax in aldi eranoe 


m aXciisrnibh* (Q ) ^ iVnsrmbU* (Q ) /q a) 

Q= (Q) “ 

is given for various quan tities Q like the wake width, the resolved turbulant kinetic 
energy density integrated in y , and the resolved turbulant kinetic energy disspation 
rate integrated in y . For these global quantities, Q is originally 0 and rapidfy reaches 
values of the order of 5% .We investigated the behaviour of Q fora local quantity, the 
m ax in urn griipoint value of the stream w ise vorticity com ponent u) x ^ uix . H ere, Q is 
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FIGURE 5. Test of de-correlation of the LES fields. The reduced maximal difference A Q is plotted 
versus time for the wake width (left, solid line), the resolved turbulent kinetic energy density 
integrated in y (left, dashed line), the resolved turbulent kinetic energy dissipation integrated 
in y (left, dotted line) and the maximum grid-point value of the .r-component of the resolved 
vorticity, u/ T , r liax (right) 

origin ally 5% and rapid V reaches values of the order of 20 to 40% . T hese m easuren cuts 
indicate a fairV rapid de-<torre3ation of the di erent m an bers of the ansan ble. 

A gain, we stress here that thism ethodology used for building the initfol conditions has 
the m a jor advantage of guaranteeing that the statisticalprcperties are initially identical 
for each realization, whife results p respited in gure 5 indicate reasonable de-oorreh tons- 
il ow ever, w hen no D N S is available, it w ill be unavoidable to develop a suitable approach 
for building in ithl conditions thatw ill satisfy the criteria of independence and equivalence 
as described in Section 2d. In that case, the constraints to be satis ed by the initi al 
conditions w dim ost probably oom e firm experin an tal data. 

5 3. T es ts of convergence 

In order to test the convergence of the EAD P results for increasing values of R , two 
types of teste were perform ed. First, the an san ble-averaged values of several relevant 
quantities in the tin e-evolving wake cm have bean com pared for various ansan ble sizes. 
Tn particular, the results for (z) the w ake w idth , (n) the turbulait kinetic energy density 
integrated in y and (ziz ) the turbulent kinetic energy dissipation rate integrated in y are 
compared for R - 4, 8, 16, and 32. As can be seen firm gure 6, only the turbulent 
kinetic energy integrated in y is a ected by the num ber of realisations. H cwa/er, the 
values obtained with 16 and 32 realisations are aln ost indistinguishable for all three 
quantities. 

Second, the in uence of the ansan ble size on the oom puted eddy viscosity was exam - 
ined . The pro le of the m ean eddy viscosity and the faction of grid points for which 
the eddy viscosity has been clpped according to the criterion (5.6) are ccm pared for 
the sam e values of R in gures 7 and 8. As seen in gure 7 , the eddy viscosity pro le 


1 





Wake width 


17 


S tatis tkal ens an ble of hrge eddy s in uhtfons 



Time 

FicuiRE fi. Convergence of the enseinbie-averageci evolutions of the wake width (top left), the 
resolved turbulent kinetic energy density integrated in y (top right) and the resolved turbulent 
kinetic energy dissipation integrated in y (bottom). Various ensemble sizes tire compared: Ft = 4 
o, R = 8 □, R = 16 o, and R = 32 A. 


depends onV weak V on the num ber of realisations forvaliesof/? between. 4 and 32, and 
the pro les are n ea rly identical for R = 16 and i? = 32 . A s expected , the fraction of grid 
points requiring clpphg of the m odel coe ciant C rapid V decrease w ith R { gure 8 ) . 
T he total fraction of clpped points iutegrated h y is less than 1% for R = 16 during the 
entire sin ulatdon . T his, can hired w ith the very an all change in m ost of the ensenble- 
averaged quantities as R is increased foam 16 to 32, supports the adoption of R - 16 
as a reason abfe ausem ble size for both m odel testhg and production LE S . B ecause this 
value of R is the sarn e as that required for the sin u&tion of hem ogeneous turbulence, it 
333 * s reasonable to hope that R= 16 provides an adequate arsem ble size for the E A D P 
in evar m ore ccm plicated geem e tries. 

T he can parison betwear various arson ble sizes is presented here onV for m odel A 
(2 A ) . T he sam e conclusions concerning the convergence of the results and the appropriate 
value of Z? are obtained when either the Sm agorinsky m odel or m odel B (25) is used. 
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Figure 7. Convergence of the eddy viscosity profile for /. = 0 (left) ;ind for t = 250 (right). 
Various ensemble sizes are compared: i? = 4o, 7? = 8o,/? = lGo, and R = 32 A. The eddy 
viscosity is normalized by the molecular viscosity. 




FIGURE 8. Profile of the fraction of grid points requiring clipping of the coefficient ( ' at t = 0 
(left) and / = 250 (right). Two ensemble sizes are compared: R — 4 solid line and R — 32 dashed 
line. 


5 .4 . C cm parison of m odeJs 

A s already m action ed , an in portant m otiyatbn for developing the E A D P is the possi- 
bility of investigating new concepts in subgrid -scale m odelling. H ere, the Itered DNS 
of M oser, R ogees, & Ew ing (1997) is can pared w ith the LE S predictions of G hosal & 
R ogers (1997) and the predictions of the m odels presented in section 51.W e have also 
added the results of a LE S without a sub grid -scale m odeL In all cases, and in agrean ant 
w ith the conclusion of the preceding section, the sin u&tions for the EAD P have bean 
perform ed w ith R - 16 . 

T he rst in portant conclusion is that the plane-averaged and ansan ble-averaged dy- 
nam ic procedures fcad to iidistinguidnable results whan they are applied w ith the sam e 
m odeL For instance, in the LES of G hosal & R ogers, the plane-averaged dynam ic proce- 
dure w as in plan anted w ith the standard Sm agorindey m odeL T heir results are identical 
to those obtained whan the Sm agorinsky m odel is used with the EAD P . In the follow ing 
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Figure 9. Evolution of the wake width (top left), the resolved turbulent kinetic eueig> density 
integrated in y (top right) and the resolved turbulent kinetic energy production integrated m y 
obtained from the filtered DNS the Smagorinsky model — , the model A □, the model B + 
and no model x . 

can panson , the Sm agorinsky case w ill refer to both the EA D P and the plane-averaged 
LES ofGhosalfc R ogers. 

T he evohtbns of {i ) the w ake w idth f (n) the resolved turbubnt kinetic energy density 
integrated in y , and (nz ) the resolved turbulent kinetic oiergy production integrated in 
y aje presented in gure 9 . The wake width is dcm hated by large-scale ow features 
and consequently is not strongly a ected by the models. Actually, the prediction of 
the LE S w ithout a subgrid -scale m odel {an under-resolved DNS) provides a reasonable 
approach ation to the value obtained by Itering the D N S data. 

T he turbufent kinetic energy density integrated in y is m ore di cult to predict using 
LES. N ot using a subgrid -scale m odel results in poor prediction of resolved energy dan- 
sity.M odelA leads to aln ost the sam e result as the Sm agorhsky m odel This is a ganeral 
feature of the dynam ic procedure that has bean noted previously (W ong & Lilly 1994, 
Cara tret al 1995b) . H ow ever, w ith in the dynam ic procedure approach, m odel A is can - 
pu tatbnalV m uch cheaper to in plan ant than the Sm agorhsky m odel and this m otivates 
the use of the scalhg (52) for the sub -grid stress instead of (41). Them odel B , fern 
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which the mean ble-averaged resolved strain aate has bean ran oved, leads to results that 
better t the DNS data in the earV stages of the sin ulation . A t later Lm es, however, 
this m odel is farther ffan the Itered DNS valies than m odelA and the Sm agorin^cy 
m ode! Inganeral the predictions of all three m odels sean comparable. 

T he evolution of the turbulent kinetic energy production also dncw s the in portant role 
of the m odels. The no-m odelLES prediction for the resolved energy production is m uch 
too high in the early stage and too lew at later tin es. A gain , m odelA leads to aim ost the 
sam e result as the Sm agorinsky m odel M odelB system atically over-p red acts the energy 
production. H aw ever, it would be rather qpecuhtive to draw any de native conclusion 
regarding which m odel (A or B ) is better ham the results presented here. 

6 . C ondnsions 

The teste presented here have dnewn that the knowledge of statistically equivalent 
resolved velocity elds m ay be useful in deriving new subgrid -scale m odels. W e have used 
the additional inform a toon available firm the di erant LES 's to propose an ansam ble- 
averaged version of the dynam ic procedure. This dynam oc procedure presents several 
advantages. 

F irst, a local version of the ansan ble-averaged dynam ic m odel is derived in the lin it of 
large ansan ble sets. T he local form ulation does not rely on any hem oganeity assum ption 
and can be adapted to any geom etry, unlike the chissical ^oiim e- (or plane- or line-) 
averaged dynam ic m odels. 

Second, the practical lin it of large an sen ble sets is closely approached for R zs 16. 
This is indicated by m any diagnostics. For example, the PDF of the m odelparam eter 
appears to be very peaked for R = 16 and its qoatal vara tons decrease drastically for 
increasing ansan ble sizes and saan to be quite m did for R - 16 . A Iso, all the m easured 
quan tities , both in horn ogaieous turbulence and in the plane wake, are ain ost identical 
for R = 16 and R = 32 . T his is, of course, a m apr an couragan ant for farther developing 
the E A D P m ethodology . T he fact that the sam e value R - 16 appears to be appropriate 
for both hem oganeous turbu lance and the plane wake suggests that this m ight be an 
adequate an sen ble size for converged results even in m ore com plex ow s. 

A Iso, the E A D P reconciles the dynam ic procedure w ith the concept of a universal 
param eter in turbu lance m odelling . T his is a very desirable property since it is can m only 
accepted that m odelparam eters dnould depend onV on the external conditions of the 
cm and not on the particular realisation that is observed . 

C on sidering the rapid developm ant of parallel ocm puters, the use of an an sen ble of 
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statistically equivalent and hdqpaidait LE S 's can be regarded as a very prom isirg tech- 
nique. T his technique can be in plen ented w ith feirly 9ti all ersan ble sizes. 0 riginal 
m odelling concepts that cannot be in plan anted in fiilV xihcm ogeneous cwsby con- 
vantbnalLES techniques are possible within the fram ework of the EADP and warrant 
further exam hatbn . M oreover, the additional ocm putational cost generated by the use 
of R sin ultaneous LES's oould be ccm pensated by the feet that statistics can be accu- 
m u b ted m uch fester w ith the arisen ble of LE S 's if one ensures that all the realisations 
are statistically independent. 

This work was conpleted prm ariV during the 1996 and 1998 summ er program s or- 
ganized by the Center for Turbulence Research, Stanford University & NASA Am es 
R esearch C enter. FruitfhlD iscussionsw ith W ill am C abot and P arviz M oh are ackncw 1- 
edged. D C . is \Chercheur Quali e du Fonds National de la Recherche Scion ti que" , 
Barium .This work has bean supported by the \ convention FR FC 2.4563 38" and by the 
NATO grant C RG 970213. 
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